Search for: All records

About 70 luminous quasars discovered atz> 6.5 are strongly biased toward the bright end, thus not providing a comprehensive view of quasar abundance beyond the cosmic dawn. We present the predicted results of the Roman/Rubin high-redshift quasar survey, yielding 3 times more, 2–4 mag deeper quasar samples, probing high-redshift quasars across a broad range of luminosities, especially faint quasars atL_bol∼ 10¹⁰L_⊙orM₁₄₅₀∼ −22, which are currently poorly explored. We include high-zquasars, galactic dwarfs, and low-zcompact galaxies with similar colors as quasar candidates. We create mock catalogs based on population models to evaluate selection completeness and efficiency. We utilize the classical color dropout method in thezandYbands to select primary quasar candidates, followed up with the Bayesian selection method to identify quasars. We show that overall selection completeness >80% and efficiency ∼10% at 6.5 <z< 9, with 180 quasars atz> 6.5, 20 atz> 7.5, and 2 atz> 8.5. The quasar yields depend sensitively on the assumed quasar luminosity shape and redshift evolution. Brown dwarf rejection through proper motion up to 50% can be made for stars brighter than 25 mag, low-zgalaxies dominate at fainter magnitude. Our results show that Roman/Rubin are able to discover a statistical sample of the earliest and faintest quasars in the Universe. The new valuable data sets are worth follow-up studies with JWST and Extremely Large Telescopes to determine the quasar luminosity function faint end slope and constraint the supermassive black holes growth in the early Universe.

 

Protoclusters, the progenitors of galaxy clusters, trace large scale structures in the early Universe and are important to our understanding of structure formation and galaxy evolution. To date, only a handful of protoclusters have been identified in the Epoch of Reionization. As one of the rarest populations in the early Universe, distant quasars that host active supermassive black holes are thought to reside in the most massive dark matter halos at that cosmic epoch and could thus potentially pinpoint some of the earliest protoclusters. In this Letter, we report the discovery of a massive protocluster around a luminous quasar atz= 6.63. This protocluster is anchored by the quasar and includes three [Cii] emitters atz∼ 6.63, 12 spectroscopically confirmed Lyαemitters (LAEs) at 6.54 <z≤ 6.64, and a large number of narrow-band-imaging selected LAE candidates at the same redshift. This structure has an overall overdensity of$\delta ={3.3}_{-0.9}^{+1.1}$within ∼35 × 74 cMpc²on the sky and an extreme overdensity ofδ> 30 in its central region (i.e.,R≲ 2 cMpc). We estimate that this protocluster will collapse into a galaxy cluster with a mass of${6.9}_{-1.4}^{+1.2}\times {10}^{15}{M}_{\odot }$at the current epoch, more massive than the most massive clusters known in the local Universe such as Coma. In the quasar vicinity, we discover a double-peaked LAE, which implies that the quasar has a UV lifetime greater than 0.8 Myrs and has already ionized its surrounding intergalactic medium.

 

Using recent empirical constraints on the dark matter halo–galaxy–supermassive black hole (SMBH) connection from z = 0–7, we infer how undermassive, typical, and overmassive SMBHs contribute to the quasar luminosity function (QLF) at z = 6. We find that beyond Lbol = 5 × 1046 erg s−1, the z = 6 QLF is dominated by SMBHs that are at least 0.3 dex above the z = 6 median M•–M* relation. The QLF is dominated by typical SMBHs (i.e. within ±0.3 dex around the M•–M* relation) at Lbol ≲ 1045 erg s−1. At z ∼ 6, the intrinsic M•–M* relation for all SMBHs is slightly steeper than the z = 0 scaling, with a similar normalization at $M_* \sim 10^{11} \, \mathrm{M}_\odot$. We also predict the M•–M* relation for z = 6 bright quasars selected by different bolometric luminosity thresholds, finding very good agreement with observations. For quasars with Lbol > 3 × 1046 (1048) erg s−1, the scaling relation is shifted upwards by ∼0.35 (1.0) dex for 1011M⊙ galaxies. To accurately measure the intrinsic M•–M* relation, it is essential to include fainter quasars with Lbol ≲ 1045 erg s−1. At high redshifts, low-luminosity quasars are thus the best targets for understanding typical formation paths for SMBHs in galaxies.

 

We present the first results from a new survey for high-redshift (z≳ 5) gravitationally lensed quasars and close quasar pairs. We carry out candidate selection based on the colors and shapes of objects in public imaging surveys, then conduct follow-up observations to confirm the nature of high-priority candidates. In this paper, we report the discoveries of J0025–0145 (z= 5.07), which we identify as an intermediately lensed quasar, and J2329–0522 (z= 4.85), which is a kiloparsec-scale close quasar pair. The Hubble Space Telescope (HST) image of J0025–0145 shows a foreground lensing galaxy located 0.″6 away from the quasar. However, J0025–0145 does not exhibit multiple lensed images of the quasar, and we identify J0025–0145 as an intermediate lensing system (a lensing system that is not multiply imaged but has a significant magnification). The spectrum of J0025–0145 implies an extreme Eddington ratio if the quasar luminosity is intrinsic, which could be explained by a large lensing magnification. The HST image of J0025–0145 also indicates a tentative detection of the quasar host galaxy in the rest-frame UV, illustrating the power of lensing magnification and distortion in studies of high-redshift quasar host galaxies. Object J2329–0522 consists of two resolved components with significantly different spectral properties and a lack of lensing galaxy detection under subarcsecond seeing. We identify it as a close quasar pair, which is the highest confirmed kiloparsec-scale quasar pair to date. We also report four lensed quasars and quasar pairs at 2 <z< 4 and discuss possible improvements to our survey strategy.

 

We present band 6 ALMA observations of a heavily obscured radio-loud (L1.4 GHz = 1025.4 W Hz−1) active galactic nucleus (AGN) candidate at zphot = 6.83 ± 0.06 found in the 1.5 deg2 COSMOS field. The ALMA data reveal detections of exceptionally strong [C ii]158 $\mu$m (z[C ii] = 6.8532) and underlying dust continuum emission from this object (COS-87259), where the [C ii] line luminosity, line width, and 158 $\mu$m continuum luminosity are comparable to those seen from z ∼ 7 sub-mm galaxies and quasar hosts. The 158 $\mu$m continuum detection suggests a total infrared luminosity of $9\times 10^{12}\, \mathrm{ L}_\odot$ with corresponding very large obscured star formation rate (1300 M⊙ yr−1) and dust mass ($2\times 10^9\, \mathrm{ M}_\odot$). The strong break seen between the VIRCam and IRAC photometry perhaps suggests that COS-87259 is an extremely massive reionization-era galaxy with $M_\ast \approx 1.7\times 10^{11}\, \mathrm{ M}_\odot$. Moreover, the MIPS, PACS, and SPIRE detections imply that this object harbours an AGN that is heavily obscured ($\tau _{_{\mathrm{9.7\,\mu m}}}=2.3$) with a bolometric luminosity of approximately $5\times 10^{13}\, \mathrm{ L}_\odot$. Such a very high AGN luminosity suggests that this object is powered by an ≈1.6 × 10$^9\, \mathrm{ M}_\odot$ black hole if accreting near the Eddington limit, and is effectively a highly obscured version of an extremely ultraviolet (UV)-luminous (M1450 ≈ −27.3) z ∼ 7 quasar. Notably, these z ∼ 7 quasars are an exceedingly rare population (∼0.001 deg−2), while COS-87259 was identified over a relatively small field. Future very wide area surveys with e.g. Roman and Euclid have the potential to identify many more extremely red yet UV-bright z ≳ 7 objects similar to COS-87259, providing richer insight into the occurrence of intense obscured star formation and supermassive black hole growth among this population.

 

We present the astrometric calibration of the Beijing–Arizona Sky Survey (BASS). The BASS astrometry was tied to the International Celestial Reference Frame via the Gaia Data Release 2 reference catalog. For effects that were stable throughout the BASS observations, including differential chromatic refraction and the low charge transfer efficiency of the CCD, we corrected for these effects at the raw image coordinates. Fourth-order polynomial intermediate longitudinal and latitudinal corrections were used to remove optical distortions. The comparison with the Gaia catalog shows that the systematic errors, depending on color or magnitude, are less than 2 milliarcseconds (mas). The position systematic error is estimated to be about −0.01 ± 0.7 mas in the region between 30° and 60° of decl. and up to −0.07 ± 0.9 mas in the region north of decl. 60°.

 

We measure the mean free path ($\lambda _{\rm mfp,H\, \small {I}}$), photoionization rate ($\langle \Gamma _{\rm H\, \small {I}} \rangle$), and neutral fraction ($\langle f_{\rm H\, \small {I}} \rangle$) of hydrogen in 12 redshift bins at 4.85 < z < 6.05 from a large sample of moderate resolution XShooter and ESI QSO absorption spectra. The fluctuations in ionizing radiation field are modelled by post-processing simulations from the Sherwood suite using our new code ‘EXtended reionization based on the Code for Ionization and Temperature Evolution’ (ex-cite). ex-cite uses efficient Octree summation for computing intergalactic medium attenuation and can generate large number of high resolution $\Gamma _{\rm H\, \small {I}}$ fluctuation models. Our simulation with ex-cite shows remarkable agreement with simulations performed with the radiative transfer code Aton and can recover the simulated parameters within 1σ uncertainty. We measure the three parameters by forward-modelling the  Lyα forest and comparing the effective optical depth ($\tau _{\rm eff, H\, \small {I}}$) distribution in simulations and observations. The final uncertainties in our measured parameters account for the uncertainties due to thermal parameters, modelling parameters, observational systematics, and cosmic variance. Our best-fitting parameters show significant evolution with redshift such that $\lambda _{\rm mfp,H\, \small {I}}$ and $\langle f_{\rm H\, \small {I}} \rangle$ decreases and increases by a factor ∼6 and ∼104, respectively from z ∼ 5 to z ∼ 6. By comparing our $\lambda _{\rm mfp,H\, \small {I}}$, $\langle \Gamma _{\rm H\, \small {I}} \rangle$ and $\langle f_{\rm H\, \small {I}} \rangle$ evolution with that in state-of-the-art Aton radiative transfer simulations and the Thesan and CoDa-III simulations, we find that our best-fitting parameter evolution is consistent with a model in which reionization completes by z ∼ 5.2. Our best-fitting model that matches the $\tau _{\rm eff, H\, \small {I}}$ distribution also reproduces the dark gap length distribution and transmission spike height distribution suggesting robustness and accuracy of our measured parameters.

 

ABSTRACT We introduce a probabilistic approach to select 6 ≤ $z$ ≤ 8 quasar candidates for spectroscopic follow-up, which is based on density estimation in the high-dimensional space inhabited by the optical and near-infrared photometry. Densities are modelled as Gaussian mixtures with principled accounting of errors using the extreme deconvolution (XD) technique, generalizing an approach successfully used to select lower redshift ($z$ ≤ 3) quasars. We train the probability density of contaminants on 1902 071 7-d flux measurements from the 1076 deg2 overlapping area from the Dark Energy Camera Legacy Survey (DECaLS) ($z$), VIKING (YJHKs), and unWISE (W1W2) imaging surveys, after requiring they dropout of DECaLS g and r, whereas the distribution of high-$z$ quasars are trained on synthetic model photometry. Extensive simulations based on these density distributions and current estimates of the quasar luminosity function indicate that this method achieves a completeness of $\ge 56{{\ \rm per\ cent}}$ and an efficiency of $\ge 5{{\ \rm per\ cent}}$ for selecting quasars at 6 < $z$ < 8 with JAB < 21.5. Among the classified sources are 8 known 6 < $z$ < 7 quasars, of which 2/8 are selected suggesting a completeness $\simeq 25{{\ \rm per\ cent}}$, whereas classifying the 6 known (JAB < 21.5) quasars at $z$ > 7 from the entire sky, we select 5/6 or a completeness of $\simeq 80{{\ \rm per\ cent}}$. The failure to select the majority of 6 < $z$ < 7 quasars arises because our quasar density model is based on an empirical quasar spectral energy distribution model that underestimates the scatter in the distribution of fluxes. This new approach to quasar selection paves the way for efficient spectroscopic follow-up of Euclid quasar candidates with ground-based telescopes and James Webb Space Telescope. 

ABSTRACT We study the demographics of z ∼ 6 broad-line quasars in the black hole (BH) mass–luminosity plane using a sample of more than 100 quasars at 5.7 < z < 6.5. These quasars have well-quantified selection functions and nearly one-third of them also have virial BH masses estimated from near-IR spectroscopy. We use forward modelling of parametrized intrinsic distributions of BH masses and Eddington ratios, and account for the sample flux limits and measurement uncertainties of the BH masses and luminosities. We find significant differences between the intrinsic and observed distributions of the quantities due to measurement uncertainties and sample flux limits. There is also marginal evidence that the virial BH masses are susceptible to a positive luminosity-dependent bias (BH mass is overestimated when luminosity is above the average), and that the mean Eddington ratio increases with BH mass. Our models provide reliable constraints on the z ∼ 6 BH mass function at $M_{\rm BH}\gt 10^{8.5}\, M_\odot$, with a median 1σ uncertainty of ∼0.5 dex in abundance. The intrinsic Eddington ratio distribution of $M_{\rm BH}\gt 10^{8.5}\, M_\odot$ quasars can be approximated by a mass-dependent Schechter model, with a broad peak around log (Lbol/LEdd) ∼ −0.9. We also find that, at 4.5 ≲ z ≲ 6, the number densities of more massive BHs tend to decline more rapidly with increasing redshift, contrary to the trend at 2.5 ≲ z ≲ 4.5 reported previously. 

Abstract We present a well-designed sample of more than 1000 type 1 quasars at 3.5 < z < 5 and derive UV quasar luminosity functions (QLFs) in this redshift range. These quasars were selected using the Sloan Digital Sky Survey (SDSS) imaging data in the Stripe 82 and overlap regions with repeat imaging observations that are about 1 mag fainter than the SDSS single-epoch data. The follow-up spectroscopic observations were conducted by the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) as one of the BOSS ancillary programs. Reaching i ∼ 21.5 mag, our sample bridges previous samples from brighter and deeper surveys. We use a 1/ V a method to derive binned QLFs at 3.6 < z < 4.0, 4.0 < z < 4.5, and 4.5 < z < 4.9 and then use a double power-law model to parameterize the QLFs. We also combine our data with literature QLFs to better constrain the QLFs across a much wider luminosity baseline. The faint- and bright-end slopes of the QLFs in this redshift range are around −1.7 and −3.7, respectively, with uncertainties from 0.2 to 0.3 to >0.5. The evolution of the QLFs from z ∼ 5 to 3.5 can be described by a pure density evolution model (∝10 kz ) with a parameter k similar to that at 5 < z < 7, suggesting a nearly uniform evolution of the quasar density at z = 3.5–7.